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Looking inside the nest: the hidden structure of the merger galaxy 
NGC 1316 (Fornax A)* 
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ABSTRACT 

We present an analysis of the circumnuclear structure of NGC 1316 using both near-infrared 
imaging and stellar kinematics. 2D decomposition of the images suggests the presence of a 
structure that resembles inner gaseous spiral arms, at about 5 to 15" from the center (ss 500 
to 1500 pc). We also find a disc-like stellar structure with radius less than 200 pc. Analy- 
sis of previously published SINFONI integral field kinematics data indicates a kinematically 
decoupled core in the same spatial scale, further evidence that indeed the nuclear stellar struc- 
ture found is a kinematically cold stellar disc. We suggest that both newly-found structural 
components are the result of a recent accretion of a companion galaxy. 

Key words: galaxies: evolution - galaxies: formation - galaxies: individual: 
(NGC 1316=Fornax A) - galaxies: interactions - galaxies: peculiar - galaxies: structure 



; 1 INTRODUCTION 

NGC 1316 (Fornax A) is a giant elliptical radio galaxy (Fig. Q] 
- top panels) wi th pronounced dust pa tches, Ha filaments, rip- 
1 pies, and loops dSchweizerlll980l . Il98lh . located in the outskirts 
of the Fornax cluster, 3? 7 away from the central giant ellipti- 
cal NGC 1399. From its high luminosity and low central veloc- 
ity dispe rsion compared to that expected from the Faber- Jackson 
relation dD'Onofrio et"aI1l 19951). NGC 1316 has been considered 
as a merger remnant teosma etallll985h with a tidal-tail system 
comprising five tails or loops of varying morph ology. A merger 
event seemed to have occurred about 3 Gyr ago dGoudfrooii et al.l 
1200 lh . while evidence for two more recent mergers have also been 
found dSchweized[l98(il ; iMackie and Fabbianoll 19981) . Cataloged 
as a SO peculiar galaxy, it has been extensively observed in a 
wide range of wavelengths, and has non-stellar nuclear activity. 
In the r adio domain, Forn ax A is one of the brightest objects in 
the sky dEkers et alj[l983h . It has prom inent giant radio lobes at 
a position angl e of ^110 ° dWad el ll96ll) consisting of polarized 
filaments dFomalont et al.l 19891) and S-shaped nuclear radio jets 
dGeldzahler and Fomalonj|l984f) . Stellar kinematics of NGC 1316 
shows that the gala xy is rotating around the minor axis at a po sition 
angle of ~ 140° dBosma et alll 19851 ; lArnaboldi et al.ll 19981) . The 
kinematics of the outer halo of NGC 1316 indicates that this early- 
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type galaxy contains as much angular mom entum as a giant spiral 
of similar luminosity ( Arnaboldi et al. 1998). 

Decomposition of NGC 1316 radial surface brightness pro- 
file has been performed in the past. The excess of brightness in 
the central region (deviations from r 1,/4 law) has been noted in 
ISchweized dl98Cl) . iD'Onofricl d200ll) found a nuclear disc from 
ground-based B-band observations, while Ide Souza et all (2004), 
also using ground-based optical data, obtained different parameters 
for this disc. 

In this paper, we revisit the central region of NGC 1316. Using 
ground-based near-infrared images and kinematics we show that 
the discovered brightness excess in the central region of the galaxy 
(within 2 arcsec) might be explained by the presence of a kinemat- 
ically decoupled core. 

For the remaining of this paper, we adopt a distance D = 
18.6 Mpc, based o n HST measurements of Cepheid variables 
dMadoreetai]|l999l) in NGC 1316, so that 1 arcsec corresponds 
to 90 pc. 



2 OBSERVATIONS AND DATA REDUCTION 

The photometric observations were made on the 7th and 8th of 
March 2001 using the near-IR imager/spectrometer SOFI (J, H, 
and Ks bands) on the 3.5 m NTT telescope in La Silla, Chile. The 
SOFI detector is a Hawaii HgCdTe 1024 x 1024 array. We used 
the SOFI large-field mode, with a pixel scale 0.29 arcsec per pixel 
and a field of view 4.9 arcmin. Conditions during the observations 
were generally clear with a typical seeing better than 0.7 arcsec. 
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Figure 1. SOFI Ks image of NGC 1316 emphasizing the outer (top left panel, a) and inner (top right panel, b) parts of the galaxy. The middle panels (c,d) 
show the residual image after subtracting the BUDDA bulge+agn model, while the bottom panels (e,f) show the residual image from the bulge+disc+agn model 
and the spatial scale is the same as in the corresponding top panels. The color coding is indicated in the right-hand bars, in magnitudes per square arcsecond. 
The square in panels c and f represents the SINFONI field of view. The inner gaseous spiral arms are seen as negative residuals (yellow and red) in panel f. 
The inner stellar disk is seen as positive residual (blue and green), interior to the gaseous spirals and in the center of the SINFONI field of view. (The very 
central negative residual is 2 pixels in size, and it is smaller than the seeing FWHM.) The two arrows in panel f show positions in the inner and outer parts 
of the gaseous spirals where the residuals are strongly negative. They coincide with remarkable features seen in the geometric profiles of Fig. [2] which are 
consistent with the presence of inner gaseous spiral arms. (North is up, East to the left.) 
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Figure 2. Results from ellipse and BUDDA fits to the SOFI Ks-bund im- 
age. The top left panel shows the surface brightness radial profiles of the 
galaxy (dots with error bars) and the BUDDA models with and without a 
disc component (blue/solid and red/dashed lines, respectively - both models 
also include an AGN component). The top right panel is the radial profile 
of ellipticity, and the corresponding profiles of position angle and the 64 
Fourier coefficient are shown at bottom left and right panels, respectively. 
The dashed and solid arrows show the position of the features seen in the 
residual image in Fig. [if. The thin magenta lines are produced from ellipse 
fits to an image of the disc+bulge model, in which the model inner spirals 
arms seen in Fig. [3] are subtracted (only the inner part of the ellipse fits is 
shown - see text for details). They reproduce well the inner features seen 
in the galaxy geometric profiles, suggesting the presence of inner gaseous 
spirals in NGC 1316. 

The integration times were 20 x 6 s for each band. The observa- 
tions were carried out using the "jitter" technique. Jittering was 
controlled by an automatic jitter template, which produces a set of 
dithered frames. Offsets are generated randomly within a box of 
40 x 40 arsec centered on the galaxy center. Since the full size 
of the galaxy observed is bigger than the field of view of SOFI 
we obtained intermediate sky frames of the same exposure time. A 
classical near-infrared data reduction procedure was adopted. Dark 
frame subtraction, flat fielding and illumination correction were ap- 
plied using CCDRED tasks in IRAlfl The final sky-subtracted im- 
ages were aligned and combined using JITTER rou tine of the ESO's 
ECLIPSE package (version 4.9.0, Devi liard|[l997h . The zero-points 
were checked by comparing the photometry of stars with measure- 
ments from 2MASS. 



3 2D IMAGE DECOMPOSITION 

We used BUDDA jde Souza et alj|2004 iGadottil feOQS) to produce 
2D image models of NGC 1316, using the SOFI ife-band image. 
The PSF is modeled using several suitable stars in the field and us- 
ing a Moffat function, which is k nown to produce better results than 
a Gaussian hruiillo et ai]|200ll) . The first model corresponds to a 
single component, described with a Sersic surface brightness radial 
profile and fixed position angle and ellipticity, plus a point source, 
convolved with the PSF, to model the AGN. The resulting fit indi- 
cates a Sersic index n — 4.2 and an effective radius of 96 arcsec. 
Since the inner half of the galaxy seems somewhat boxy, we left the 

1 Image Reduction and Analysis Facility (IRAF) is distributed by NOAO, 
which is operated by AURA, Inc., under contract to the NSF. 



30" 




SPIRAL MODEL 



Figure 3. Model inner spiral arms built with IRAF and used to simulate the 
effects of inner gaseous spiral arms on radial profiles of geometric parame- 
ters from ellipse fits to the disc+bulge model of NGC 1316 (see the magenta 
lines in Fig. [2}- The color code is the same as in Fig. [TJ. As in Fig. [J, the 
inner square shows the SINFONI field of view. (North is up, East to the 
left.) 



boxiness parameter c free, but that resulted in c = 2, correspond- 
ing to a perfect ellipse (i.e. neither boxy or disky). We obtained a 
residual image by subtracting this model from the original galaxy 
image, and verified that a significant inner structure is seen (middle 
panels of Fig. Q}. We produced test fits with an offset centroid to 
check whether this feature could be a result of an ill-defined cen- 
troid, but we find that the uncertainty in the centroid position is of 
1 pixel only. In any case, a bad centroid position would produce a 
region of negative residual alongside another of positive residual, a 
pattern which is different from what is seen in Fig. QJi. This panel 
shows an oscillation between negative and positive residuals. This 
also cannot be the product of a mismatch between the PSF in the 
image and how we modeled the PSF, as the residuals are an order 
of magnitude larger than the PSF. This motivated us to produce a 
second model, which includes a disc component following an ex- 
ponential surface brightness radial profile. The best solution model 
in this case has a disc component with a scale length of 58 arcsec, 
and a bulge component with Sersic index n = 2.9 and effective 
radius of 17 arcsec. The bulge component is now slightly boxy, 
with c = 2.2, and the bulge/total ratio is 0.37. This suggests that 
by adding another component the code is able to retrieve better the 
boxy aspect of the bulge. This model produces a much better fit 
to the galaxy image, as can be seen in Fig. [2] which shows results 
from ellipse fits to the galaxy image and the models, using IRAF. 
The reduced \ 2 g° es down from 5.5 to 3.8 with a total of 12 pa- 
rameters fitted, instead of 8 (4 new parameters to describe the disc). 
Further, the decomposition with only the bulge model results in a 
residual image with positive residuals amounting to ~ 16% of the 
total galaxy light. 

The residual image from the first BUDDA model (Fig.[T]l) and 
the ellipse fits (Fig. [2} give a strong argument for the presence of 
an extra inner component in NGC 1316. The difficulty in interpret- 
ing such substructure is that it does not resemble morphologically a 
well-know stellar structure. This might be a result of the inadequate 
fit to the galaxy image using only one component. An extra compo- 
nent in any model usually leads to a better fit, even if such compo- 
nent in the model does not have a straightforward correspondence 
to the real structural component in the galaxy. In other words, the 
disc in the second BUDDA model might be truly fitting just an extra 
component (or components) which is not necessarily a usual disc 
as in disc galaxies. Thus, a conservative interpretation is that the 
extra disc component in the second model corresponds to substruc- 
tures which result in the morphologically perturbed appearance of 
NGC 1316 (which are well described by an exponential profile), 
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Figure 4. The SINFONI kinematic data. Velocity field of stars (left) and line-of-sight velocity dispersion map (middle) with continuum near 2.lfim contours 
superimposed. The right panel shows the cuts through the velocity and velocity dispersion maps along the major axis (PA = 53°). The systemic velocity is 
1785 km s" 1 . 



and not to a proper stellar disk. This in turn is consistent with a 
recent merger event. 

Interestingly, the residual image from the second BUDDA 
model (Fig. QJ) reveals a negative residual (i.e. a region in which 
the model is brighter than the galaxy) of less than about 0.2 mag, 
which resembles inner spiral arms. It is tempting to interpret this 
feature as inner gaseous spirals, which could be produced by the 
accretion of a companion galaxy. To verify this interpretation, we 
produced a new synthetic image as follows. We created a model of 
spiral arms resembling those seen in the negative residual of Fig. 
[TJ, and subtract it from the BUDDA model image in our second 
fit. These model inner spiral arms are shown in Fig. [3] They were 
built using ARTDATA tasks in IRAF, manually edited with the task 
IMEDIT, and convolved with a Gaussian, with the same FWHM as 
the PSF in our SOFI its-band images, using the task GAUSS. Their 
surface brightness levels have Poisson fluctuations added, decline 
roughly exponentially from the centre, and are chosen as to mimic 
the residual levels seen in Fig. QJ, i.e. below 0.2 mag per square 
arcsecond. We can thus mimic the effects of gaseous spirals atten- 
uating the light emission from the central region of NGC 1316. We 
produced ellipse fits to this new synthetic image, and the results are 
also shown in Fig. [2] (magenta lines). The ellipse fits to our new 
synthetic image reproduce well the oscillations in ellipticity, posi- 
tion angle and 64 seen in the inner part of NGC 1316. This supports 
our interpretation, suggesting the presence of inner gaseous spiral 
arms in the galaxy. It is important to note that interior to these spi- 
rals, there is a region of positive residuals of about 2", which can 
be an inner stellar disc. This will be discussed further below. 

A similar p hotometric analysis has been recently done by 
lLanz et using Spitzer archival data. They performed a 2D 

fit of the 3.6 micron image, using a bulge-only model, and find very 
similar residuals, as compared to our first fit (see their figure lc). 
However, they did not produce a fit with a disc component. They 
find structural parameters somewhat different than those from our 
first fit. They find a Sersic index of 6.1 ± 0.1, while we find a value 
of 4.2±0.9. They find an effective radius of 146±1 arcsec, whereas 
we find 96 ± 7 arcsec. It can be that a deeper Spitzer image results 
in a lager effective radius. Nevertheless, a profile with Sersic index 
6 does not differ terribly fr om one with an index of 4 (see figure 1 
in [Graham & Driver 2005). They do not show how their modeled 



profile matches the observed one. The residuals we find in Fig.QJ 
only appear when doing the second fit, including the disc. We note 
that the pixel scale and the PSF FWHM are significantly better in 
our SOFI images, as compared to the Spitzer 3.6 micron images. 
Our pixel scale is 0.29 arcsec/pix and our PSF FHWM is 0.7". The 
Spitzer image has 1.2 arcsec/pix and a PSF FWHM of 1.7". 



4 SINFONI KINEMATICS 

Information on stellar dynamics is crucial to understand the origin 
of the stellar disc-like structure suggested to be present, in the inner 
2", in the previous section. If the observed NIR brightness excess 
is really associated with the excess of stellar density we can expect 
that this region has also decoupled features in NIR kinematics. 

The most detailed study of the NGC 1316 stellar kinemat- 
ics, base d on VLT/SINFONI integral-field data, was recently pub- 
lished bv lNowaketal] (12008). Their investigation was focused on 
the very center of the galaxy, where gravitational forces of the 
supermassive black hole are important. Nina Nowak kindly pro- 
vided us with reduced SINFONI spectra taken with the largest 
field-of-view (mentioned as "250 mas data" in their original pa- 
per). The spectral data were presented in a data cube, where each 
0.25 arcsec spaxel in the 8x8 arcsec field has a galaxy spectrum 
in the range 1.95 — 2.45 pan. The stellar kinematics information 
were extracted using IDL-based software and algorithms devel- 
oped to study stellar kinemat ics with the integral-field spectrograph 
MPFS jMoiseev etai]|2004l) . The velocity field and velocity dis- 
persion map were created by cross-correlation with the spectrum 
of a K5/M0III template star (HD 181 109) in the wavelength range 
2.260 — 2.34 3/jm, which contain s high contrast 12 CO absorption 
bands. Unlike iNowak et all d2008l) . we did not use a radial or az- 
imuthal binning of the original data, in order to keep the original 
spatial sampling in the data cube. 

The maps of line-of-sight velocities (v) and velocity disper- 
sion (<r) of stars derived from the SINFONI data are displayed in 
Fig. [4] These maps reveal two kinematic features that are differ- 
ent from the simple picture expected for a rotating oblate spheroid. 
The first one is an absence of clear rotation in the central region 
(r < 1.5 — 2 arcsec). It seems as a central "plateau" of a constant 
color in the velocity field. The velocity field cross-section along the 
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galaxy major axis (see right panel in Fig.[4]l shows that the circum- 
nuclear velocity distribution is complex - a hint of low-amplitude 
counter-rotation appears on r < 1 arcsec. The second surprising 
kinematic feature is a slight central depression of the velocity dis- 
persion on r < 2 arcsec. This feature has a low contrast, but it 
is visible as a central "blue spot" on the velocity dispersion map 
and al so as a minimum i n the a radial distribution along the major 
axis. lNowak et alj J2008h discussed in detail the possible origin of 
the "cr-drop" feature using their high spatial resolution "25 mas" 
SINFONI data set. They have considered different effects (data re- 
duction artifacts, stellar population, AGN emission, etc.) and con- 
cluded that the cr-drop observed in the high resolution data could 
be caused by "either the AGN continuum emission or a cold stellar 
subsystem or both". However, they showed that an AGN distorts 
the stellar kinematics only in very center of the galaxy (r $C 0.06 
arcsec). Therefore, we suggest that the cr-drop on a larger radial 
scale is caused by a dynamically cold stellar population. We note 
that both kinematic features (an abse nce of centra l rotati on and a 
cr-drop) were already mentioned by Now ak et alj ( l2008h and our 
kinematic maps show a good agreement with their figure 6. How- 
ever, these features are more clearly visible in our un-binned maps. 
We emphasize that the spatial scale of the decoupled features in the 
stellar kinematics agrees very well with the value expected from the 
NIR brightness excess found above (r ^ 2 arcsec). Therefore, the 
connection between circumnuclear morphological and kinematical 
features is evident. 



5 DISCUSSION AND CONCLUSIONS 

A multicomponent structure in the central r egions of early- 
type galaxies is a wel l-known phenomenon teest et all l200ll ; 
lErwin and S parke 2003). Moreover, the presence of dust in the 
vicinity of central regions of post-merger galaxies (NGC 1316 is 
one of the best examples) makes a study of those regions even more 
complicated. 

A detailed analysis of the NIR structure of NGC 1316 (Figs. 
[T]and[2]l shows a negative residual suggestive of inner gaseous spi- 
ral arms from a radius of about 5" t o 15". Numerical simu lations 
of galaxy mergers, such as those by|di Matte o et al. I J2007t) . often 
show the formation of a gaseous disc with gaseous spiral arms and 
bars. We note a similar structure in the center of other giant ellip- 
tical radio galaxy, NGC 5128 (C entaurus A), observed in the mid- 
infrared bv lMirabel et all dl999i) . 

We thus put forward the interpretation that the structure we 
find near the center of NGC 1316 is indeed a gaseous component, 
most likely formed from in-falling material from a recent accre- 
tion event, which atte nuates the backgroun d stellar light from the 
galaxy. We note that lHorellou et al. (2001) detected CO, presum- 
ably associated to dust, near the center of the galaxy, although at a 
somewhat larger spatial scale than that of the structure we are inves- 
tigating here. Their CO map shows a structure which does not seem 
correlated with the negative residuals we find, but such a compari- 
son is not straightforward, as their spatial resolution (22") is rather 
coarse. iHorellouetallfeOOll) also advocate that the C O they found 
resulte d from a recent accretion. It is worth noting that lNowak et al.1 
( 2008) found a weak signature of molecular Hydrogen in a region 
northeast of the nucleus, a result which is corroborated by our own 
analysis of the SINFONI spectra. Another plausible interpretation 
is that the structure we see is caused by the galaxy AGN jet. Al- 
though we cannot rule out this possibility, we note that no emission 



lines typical of AGN are found in the SINFONI spectra, which is 
in principle inconsistent with it. 

Interior to these spirals ($: 2"), we find another structure, re- 
sembling a nuclear stellar disc (inside the square representing the 
SINFONI field of view in Fig. QJ). The stellar kinematics (Fig. 
suggests the presence of a kinematically decoupled core on the 
same spatial scale. This suggests that indeed the nuclear stellar 
structure we find in the structural analysis is a kinematically cold 
stellar disc. 

This nuclear disc could also be formed via a recent merging 
event. iD'Onofrio et al. I d 19951) noted that the central velocity dis- 
persion in NGC 1316 is lower than the one expected for a galaxy 
of this luminosity from the Faber-Iackson relation (260 versus 400 
km/s). It could be possibly explained by the presence of a dynam- 
ically cold d i sc-like component. Further evidence is presented by 
Nowa k et al. I j2008L their figure 8), which found indications of a 
younger stellar component in the very center (^ 0.5") of NGC 
1316. 
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